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HEAT TRANSFER - APPLICATION TO COMBUSTION -CHAMBER 
COOLING PASSAGES 

By Merwin Sibulkin and William K. Koffel 


SUMMARY 

A method for calculating the pressure drop of a high-speed 
compressible fluid in a constant-area duct under the simultaneous 
action of friction and heat transfer is developed. The method is 
based on the assumption of an exponential longitudinal temperature 
distribution. It is shown that the temperature distributions found 
in combustion-chamber cooling passages can be approximated by the 
assumed temperature distribution, and that a working chart based 
on this method simplifies the calculation of pressure drops across 
these passages. An illustrative example is included. 


INTRODUCTION 

An investigation of the cooling of ram-jet and tail -pipe -burner 
combustion chambers by means of air flowing through an annular cool- 
ing passage is being conducted at the NACA Lewis laboratory. The 
design of a combustion-chamber cooling system includes the calculation 
of the mass flow of cooling air necessary to maintain permissible 
wall temperatures and the determination of the pressure drop across 
the cooling passage required to obtain the necessary mass flow of 
cooling air. 

The pressure drop in a cooling passage results from the 
simultaneous action of friction and heat transfer. For low rates 
of heat transfer and a Mach number below 0.5, the pressure drop 
may be calculated by the simplified methods of references 1 to 3; 
these methods result in appreciable errors at higher Mach numbers. 
Accurate calculation of the pressure drop requires solution of a 
basic differential equation describing the variation of pressure 
of a compressible fluid under the simultaneous action of friction 
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and heat transfer. Unfortunately, this equation is not generally 
amenable to formal integration. Tables and charts that reduce the 
labor of numerical solution are given in references 4 and 5 and a 
special solution (reference 6) has been developed for the case of 
constant wall temperature. 

For one -dimensional flow, the pressure at a point is uniquely 
defined when the corresponding values of mass flow, temperature, 

Mach number, and flow area are known. The pressure change between 
two points along a constant-area duct can be determined if the temper- 
atures and Mach numbers at the two points are known. The differential 
equation relating Mach number variations and temperature distribution 
can be integrated in special cases. One of these cases occurs where 
the longitudinal temperature distribution T is given by T = ce 11 * 
(where c and n are constants and x is the distance along the 
passage). It is shown herein that the longitudinal distribution of 
air temperature in the cooling passages of ram-jet and tail-pipe- 
burner combustion chambers can be closely approximated by the use of 
the foregoing equation in a series of steps, each step having dif- 
ferent values for the constants c and n. This analysis develops 
a working chart that simplifies the determination of the Mach number 
change in a combustion-chamber cooling passage and, consequently, 
simplifies the calculation of pressure drops required across these 
passages to obtain the necessary mass flow of cooling air. The method 
developed applies to the flow of air through straight ducts having any 
constant shape, cross-sectional area, and roughness. 

The solution is presented in the form of a working chart con- 
structed for a ratio of specific heats equal to 1.40. The variables 
cover ranges sufficiently large to include almost all applications 
of engineering interest. An example illustrating the use of the 
chart for a typical problem is presented. 


SYMBOLS 


A 


B 


The following symbols are used in this report : 
flow area, square feet 
7+1 

M 2 (l+ tW +1+2rK 


1+7(1+K)M 2 


2(7+l+7K) 

7+1+27K 
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B* value of B when M = 1.0 

Db hydraulic diameter, — — — llSOL ^. r ?- a . feet 

wetted perimeter 

F friction drag, pounds 

f friction factor as defined By, — = 4f — 

A 2 D h 

f effective over-all value of friction factor f 

G mass velocity of fluid, pounds per second per square foot 

g acceleration of gravity, feet per second per second 

K . _£^h_ 

l°g e T 2 /Tl 

L* distance from any given point to point where Mach number would 
theoretically equal 1.0 

l length of passage, feet 

M Mach number 

P total pressure, pounds per square foot absolute 
p static pressure, pounds per square foot absolute 
E gas constant, foot-pounds per pound °R 
T total temperature, °E 

t static temperature, °R 

V fluid velocity, feet per second 
x distance along passage, feet 

7 ratio of specific heats 

p mass density, slugs per cubic foot 
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Subscripts : 

1 inlet 

2 outlet 


METHOD OF ANALYSIS 

For one -dimensional flow, the pressure at any point in a 
fluid is uniquely defined by the mass velocity, the temperature, 
and the Mach number at that point. The following development 
deals with the variation of Mach number along a duct. If the 
variation of Mach number is known, the corresponding pressure 
distribution is readily determined from the continuity equation 
for the one -dimensional flow of a perfect fluid as expressed by 


P 



( 1 ) 


or in terms of total temperature and total pressure 



(la) 


Equation (la) is graphically presented in figure 1. 


One form of the differential equation for the variation of 
Mach number of a compressible fluid under the simultaneous action 
of friction and heat transfer for constant area and specific heat 
is obtained from table 2 of reference 7. 



1+7M 2 





1-M 


2 



7M 2 


( ltZ r M2 ) 

1-M 2 


4fdx 

Dh 


( 2 ) 


The change in Mach number is dependent on the temperature dis- 
tribution along the duct. One of the cases for which equation (2) 
is integrable occurs when the longitudinal temperature distribution 
is given by 
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JL 

T 1 


4f x 
I> h K 


(3) 


where 


K is an arbitrary constant (fig. 2). 


Then "by differentiation 


4fcbc dT 
D^K = T 


(3a) 


If in equation (2 ) f equals f , the substitution of 
equation (3a) into equation (2) gives 


dM 2 \ 1 + it m2 . 


1 + tm 2 (1 + k) 


M 


1-M 


'J dT 
T 


Integration of equation (4) (appendix A) gives 


(4) 


/ \ ~ - 7 - + - 1 - ■ 

2 l y —1 21 7+1+27K 

2(y+l+yK) 

, ,2 7 +l+ 2 yK 

t 2 [l+y(l+K)M 2 ^J 


7+1 " Bq 


1 2 / y i 2 \ ■ 

M-l ( 1 + Ili ) 7+1+27K 


2 ( 7 +I+ 7 K) 
T] y+l+2yK 
l+7(l+K)M 1 2 

When x = l, T = Tg; then from equation (3) 


(5) 


= log -£■ = log, 
KD h e 


e Bi 


(5a) 


and 


4fL* 

KD h 


B* 


(5b) 
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where L* is the distance from any given point to the point where 
the Mach number would theoretically equal 1.0. Equation (5b) is 
graphically presented in the chart (fig. 3) for various values of K 
and for 7=1.40. (A 16 by 16 inch working chart of this figure is 
enclosed.) When K=±°°, by equation (3) Tg/T-^ = 1, which reduces 

the problem to the case of isothermal flow with friction. 

The use of the choking-length parameter 4fL*/Dh is illustrated 
in the following example: 



Given are = 0.200, K = 1.0, Z = 2.0 feet, and 4f/D^ = 0.5. 
From figure 3, 


4fL* x 


1.5 


L* x =3.0 feet 


By subtraction 


and 


L *2 = L *i " 1 “ 3.0 - 2 *° = 1»° foot 


4fL*„ 


= 0.5 
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Then, from figure 3, 


^ = 0.398 

For 7 = 1.350, the values of 4fL*/D^ are lower with a 

maximum difference of 5 percent; for 7 = 1.30, the maximum dif- 
ference would he 10 percent. 

If the Mach number along a duct is constant, according to its 
definition L* becomes infinite; if the Mach number along a duct 
decreases, its value never reaches 1.0 and L* has no positive 
value. Consequently, the results of this analysis can be applied 
only to flow systems in which the Mach number increases. 

For any initial Mach number, there is a critical value of K 
for which the Mach number along the duct is constant. This value 
can be determined by setting equation (4) equal to zero. The 
critical value of K is therefore 


K = - 1+ - 7M . 2 ( 6 ) 

7M 2 

Consequently, for a given Mach number, if the value of the friction 
factor is known, the minimum value of Tg/T^ for which the Mach 

number increases, and therefore for which the chart can be applied, 
can be determined from equations (6) and (3). 


APPLICATION TO COMBUSTION -CHAMBER COOLING PASSAGES 

An annular cooling passage formed by a concentric liner inside 
a tail -pipe burner or ram- jet combustion chamber is shown in figure 4. 
If the rise in combustion-gas total temperature is approximately 
linear, the longitudinal temperature distribution of the cooling air 
may be similar to the distributions shown in figure 5. If the initial 
cooling-air and combustion-gas temperatures are equal and no heat 
losses occur through the outer wall of the cooling passage, the 
temperature distributions will be similar to those shown for case A 
(fig. 5(a)). Such a case would be an insulated tail -pipe burner in 
which a portion of the turbine -discharge gases are used as the coolant. 
If the initial cooling-air and combustion-gas temperatures are equal, 
but the external heat losses are large, the temperature distributions 
will be similar to those shown for case B (fig. 5(b)). This case 
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would correspond to an uninsulated tail -pipe burner in which the 
outside wall of the cooling annulus is cooled by ambient air. Any 
one of the temperature distributions of case A in figure 5(a) can 
be approximated (case A, fig. 6(a)) by equation (3), where K is 
determined by letting T = Tg and x = l. The temperature distri- 
butions for case A in figure 5(a), however, have an initial slope of 
zero because the initial cooling-air and combustion-gas temperatures 
were equal. If, at the cooling-passage entrance, the cooling air 
has a lower temperature than the combustion gas, the cooling-air 
temperature distribution will have a positive initial slope. Inas- 
much as the assumed temperature distribution (given by equation (3)) 
has a positive initial slope, it is noted that for two cases (zero 
and positive initial slopes) with the same entrance Mach number and 
over-all temperature ratio, the case with a positive initial slope 
can be more closely approximated in one step. 

Any of the distributions for case A in figure 5(a) can be 
matched more closely by a two-step approximation if the temperature 
at some intermediate station (preferably near the midpoint) is 
known. In this case a separate value of K is determined for each 
step, as illustrated in appendix B. Similarly, the temperature dis- 
tributions for case B in figure 5(b) can be approached in two steps 
(fig. 6(b)). In this case an intermediate value of x near the 
minimum cooling -gas temperature would be preferable; the value of K 
for the first step would then be negative. 

The variation of Mach number (fig. 7) along the cooling passage 
for both original and approximate temperature distributions shown 
in figure 6 were calculated by a step-by-step integration of 
equation (2) using the tables of reference 4. A sufficient number 
of steps were taken to insure that any increase in the number of 
steps would not affect the result. In case A (fig. 7(a)), the 
important effect of small changes in the initial Mach number is shown. 
When the calculations show that sonic velocity is attained before the 
outlet of the cooling passage is reached, the chosen value of initial 
Mach number is greater than is physically possible, 

Case A (fig. 7(a)) shows that the Mach number increases from 
0.5 to 1.0 in a very short distance. In spite of the fact that this 
analysis is exact for the temperature distribution of equation (3), 
and can be made to closely approximate temperature distributions 
similar to those in figure_5, in a physical case small inaccuracies 
in the values chosen for f or may introduce appreciable errors 
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in the region of rapidly increasing Mach number. For example, if 
the length of the duct in figure 7(a) were 2.95 feet (fig. 8), the 
outlet Mach number based on the original temperature profile would 
be 1.0; if the value of f/D^ was then reduced by 10 percent, the 

outlet Mach number would be reduced to 0.75. At the point where 
the original value of f /D^ gives M = 0.50, however, a 10-percent 

reduction in the value of f/D^ would only reduce the Mach number 

to 0.46. These differences in outlet Mach number would give pro- 
portionate discrepancies in outlet static pressure and smaller dis- 
crepancies in outlet total pressure. Due to the uncertainty in 
choosing a friction factor or in determining the exact cross- 
sectional area for any physical case, values of outlet static 
pressure based on an outlet Mach number between 1.0 and approxi- 
mately 0.5 may therefore be subject to appreciable error. Equations 
for estimating the value of the friction factor are given in 
reference 8. The effect of the ratio of the inner to the outer 
diameter on the friction factor of an annulus is discussed in 
reference 9. 

The values of outlet Mach number obtained by the step-by- 
step integration of equation (2) along a variety of temperature - 
distribution curves similar to those of figure 5 were compared 
with the values of outlet Mach number obtained by the use of 
figure 3 (in a manner similar to that shown in appendix B), with 
the following results; 

(1) When neither method indicated choking, the outlet Mach 
number based on a one -step approximation was 0 to 25 percent greater 
than the value obtained by step-by-step integration, with the greater 
differences occurring at the higher values of outlet Mach number. 

When both methods caused choking, the theoretical, distance to choke 
based on a one -step approximation was 7 to 25 percent less than the 
distance obtained by step-by-step integration. Inasmuch as the 
temperature along these one -step-approximate curves was appreciably 
higher than the corresponding temperature on the original curves, 
this result was anticipated. 

(2) For a two-step approximation, the outlet Mach numbers 
(unchoked by either method) were within 2 percent; the difference 
in theoretical distance to cause choking (choked by either method) 
was within 2 percent. 

(3) For any value of outlet Mach number, the precentage change 
in outlet static or total pressure is equal to or less than the Mach 
number discrepancies mentioned in results 1 and 2, but in the 
opposite sense . 
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CONCLUDING REMARKS 

A method is developed for calculating the pressure drop of a 
compressible fluid under the simultaneous action of friction and 
heat transfer, and is presented in the form of a working chart. 

It is shown that the temperature distribution assumed in the 
development of the working chart can be made to closely approximate, 
in one or two steps, the temperature distributions found in the 
cooling passages of ram-jet and tail -pipe -burner combustion chambers. 
When two steps were used, the outlet pressure was within 2 percent 
of the values that would be obtained by a step-by-step integration 
along the actual temperature -distribution curve. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 17, 1949. 
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APPENDIX A 

INTEGRATION -OF -FLOW EQUATION 


One form of the differential equation for simultaneous friction 
and heat transfer (reference 7) for constant area and specific heat 
is (equation (2) of text) 



l+7M 2 (l + M 2 ) 


l-M" 


dT 

T 


7M 2 (l + M 2 


1-M £ 


4fdx 

D h 


(Al) 


If the fluid temperature distribution is given by 


dT = 4fdx 
T “ D h K 


(A2) 


substitution of equation (A2) in equation (Al) and factoring gives 


dT 

T 


(1-M 2 ) dM 2 



+ 


7-1 

2 



1+7(1+K)M 2 


(A3) 


Let a 


21li and b 
2 


7(l+K). Equation (A3) can then be written as 


dT = (l-M 2 )dM 2 

T M 2 (l+aM 2 )(l+bM 2 ) 


(A3a) 


or by separating into partial fractions, equation (A3a) becomes 


dT (l-M 2 )dM 2 a 2 (l-M^dM 2 b 2 (l-M 2 )dNf 
T " M 2 (1 + aM 2 ) + a - b (1+bM 2 ) 


(A4) 


CO 


Integration of equation (A4) results in 

'2 


^2 

r 2 

dT 

dM 2 _ 

T = 

1 J 

1 M 2 o 


A1 2. a 
dM + r — 
b-a 


V 


dM 


1 1+aM ^ 


M 2 dM 2 1 . b 2 


1+aM 


,2 r a-b 


v 


dM^ 


. 1+bM 
1 ^ 


n 2 \ 

M 2 dM 2 
1+bM 2 


(A5) 


log e T 


Jl 



2 

2 

r 

M 2 







1 

i 1 


L \ _ 1 


1+aM 2 - log e (l+aM^) 


log e (l+bM 2 ) - ^ 1+m2 - log e (l+bM 2 ) 


(A6) 



2 " 

2 

2 

2 

2 

2 ~ 

log e T 

= log M 2 
e 

1 J 

-* 2 

1 

+ S±i log (1+aM 2 ) 
^ b-a e 

1+aM 2 

b-a 

1 J 

- 1°6,(W) 

1 

L 1+bM 2 
b-a 

1 J 


(A7) 


However, 


^2 _ 1+aM^ + 1+bM^ _ 0 _ n 

b-a b-a b-a 


• * 
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Then 


log e T 


a+1 

- log M 2 (1+aM 2 )^ ~ a 


(1+bM 2 ) 


h+1 

h-a 


By substituting the limits. 


M 2 (l+aM 2 ) 


io s e = l08e 



b+1 


(l + bM 1 2 ) l '- a 


# 


Taking the antilog of equation (A9) and substituting for 
b gives 






7+1 

7+1+27K 


[l+7(l+K)M 2 


2 (7+I+7K) 

7+1+27K 




7+1 


2 \ 7+1+27K 

+ Z f M J 


1+7(1+K)M 1 2 


2 (7+1+7K) 
7+1+27K 
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(A8) 


(A9) 


and 


(A10) 
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Changing the fourth term on the right side of equation (A5) to 


+ 



P2 \ 

dM 2 M 2 dM 2 ) 

-1-bM 2 -1-bM 2 / 


and integrating in a similar manner gives a form of equation (A10) 
that can he evaluated for negative values of K. 
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APPENDIX B 
SAMPLE CALCULATION 

The pressure drop in an annular cooling passage will be cal- 
culated by the use of the large working chart in two steps. Given 
conditions : 


r 

I 


T 


1 

I 


Station * 1 1 
x = 0 
T = 500 


c d 
2.50 
603 


5.00 feet 


980° E 


Mass velocity, G, lb/sec/sq ft 14.05 

Inlet total pressure, P^, lb/sq ft absolute 1458 

Hydraulic diameter, D^, ft . . 0.0417 

Average friction factor, f. . . . 0.00756 


Inlet Mach number . - The inlet Mach number is obtained from the 
weight -flow parameter and figure 1. 


G \l^l 14.05 \J 5 OO 
T 1 = 1458 


0.2155 


Then from figure 1, for 7 = 1.40, M^ = 0.242. 

Intermediate Mach number . - The Mach number at station c is 
determined from the working chart by the following steps : 



4 X 0.00756 


0.725 


0.0417 
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4fl 


1-c 


K- 


•1-c 


. T o 
l08 el7 


0.725 x 2 .50 _ g 0y 
, 603 

log e 500 


Entering the working chart at M = 0.242 and K = 9.67 gives 


4fL* 1 


= 5.20 


which is solved for L*-],: 


L*, _ _5i20 = 7>1? feet 
1 0.725 


L * c = L* x - l 1 _ c = 7.17 - 2.50 = 4.67 feet 


Multiplication of L* by — gives 


4fL* 


0.725 X 4.67 = 3.38 


4fL* 

Reentering the working chart at - = 3.38 and K = 9.67, 

D h 

M c = 0.298 = M d . 

Outlet Mach number . - The outlet Mach number is found in the 
manner as the intermediate Mach number. 


4FZ 


d-2 


h 0.725 X 2.50 , „ 

^-2 tT ■ : 980 * 3 - 75 

loge t~ 1<>8e 605 


c 

t 

c 
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Enter the working chart at M = 0.298 and K = 3.73 and read 


4fL* d 


2.18 


Solve for L*^ 


L*, = .. 2 - 18 = 3.01 feet 
d 0.725 


and 


L *2 = L *d " Z d-2 = 3 * 01 “ 2 - 50 = °* 51 feet 


or 


4fL*. 


2 _ 


= 0.725 X 0.S1 


0.37 


Reentering the chart at 


4fL* 


''h 


0.37 and K = 3.73, M 2 = 0.567. 


Pressure drop. - The total pressure at station 2 is obtained 
from the weight-flow parameter. The weight-flow parameter at station 2 
is found from figure 1 for 7 = 1.40 and M = 0.567, 



0.431 


Solving for P 2 gives 

P = 14 •0. 5 , . 1020 lb/sq ft 

2 0.431 


The loss in total pressure across the cooling passage is 


Pi - P 2 = 1458 - 1020 = 438 lb/sq ft 


or if the loss in static pressure is desired 
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By use of reference 10, 

Pl/Pl = 0.9601 
P 2 /P 2 = 0.8040 

= 1458 X 0.9601 = 1400 lb/sq ft 

p 2 = 1020 x 0.8040 = 820 lb/sq ft 
and the loss in static pressure is 

Pl - P 2 = 580 lh/sq ft 
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Figure 1. 


Variation of weight-flow parameter 
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— - — with Mach number M. 
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Figure 2. - Temperature profile assumed in integration of differ- 
ential equation for simultaneous friction and heat transfer. 



Choking-length parameter 
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Figure 3* - Chart for determination of Mach-number variation 
in compressible fluid flowing through constant-area duct 
with simultaneous friction and heat transfer* Ratio of 
specific heats y , 1*40. (A 16- by 16-in* working chart 
of this figure is enclosed* ) 
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X = 0 



Figure 4. - Schematic diagram of combustion chamber cooled by air flowing through annular 

passage. 
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(a) Case A (no external heat losses)# 



(b) Case B (large external heat losses)# 


Figure 5. - Longitudinal distribution of cooling-air temperature in 
annular passage surrounding combustion chamber# 
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Distance x along cooling passage , ft 


(b) Case B. 

Figure 6. - Comparison of longitudinal cooling-air temperature distri- 
bution from figure 5 with distribution obtained by use of equation 
(3). 



Mach number, M Mach number 
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Distance x along cooling passage, ft 
(b) Case B. 

Figure 7. - Comparison of cooling-air Mach number distribution 
based on original and approximate temperature distributions 
of figure 6. 
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0 12 3 

Distance x along cooling 
passage, ft 


Figure 8. - Effect of on Mach number profile. 
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M=I.O 


flow area, sq ft 

hydraulic diameter ^t/ed peH^eter > ft 
friction drag, lb 

friction factor as defined in equation 

f ■* £ *, 

effective over-all value of friction factor f 

4?/ 

Oh 




theoretical distance from initial Mach number to M»I.O 

length of passage, ft 
Mach number 
total temperature, °R 
velocity, ft/sec 
mass density, slugs^u ft 



1.00 


Mach number v M 

Figure3,“ Chart for determination of Mach-number variation incompressible fiuid flowing through constant-area duct with simultaneous friction and heat transfer. Ratio of specific heats, 1.40. 
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